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Weak lensing distortion of the background cosmic microwave background (CMB) temperature and 
polarization patterns by the foreground density fluctuations is well studied in the literature. We 
discuss the gravitational lensing modification to CMB anisotropics and polarization by a stochastic 
background of primordial gravitational waves between us and the last scattering surface. While 
density fluctuations perturb CMB photons via gradient-type deflections only, foreground gravita- 
tional waves distort CMB anisotropics via both gradient- and curl-type displacements. The latter 
is a rotation of background images, while the former is related to the lensing convergence. For a 
primordial background of inflationary gravitational waves, with an amplitude corresponding to a 
tensor-to-scalar ratio below the current upper limit of ~ 0.3, the resulting modifications to the an- 
gular power spectra of CMB temperature anisotropy and polarization are below the cosmic variance 
limit. At tens of arcminute angular scales and below, these corrections, however, are above the 
level at which systematics must be controlled in all-sky anisotropy and polarization maps with no 
instrumental noise and other secondary and foreground signals. 

PACS numbers: 98.80.Es,95.85.Nv,98.35.Ce,98.70.Vc 



I. INTRODUCTION 

The weak lensing of cosmic microwave background 
(CMB) anisotropics and polarization by intervening mass 
fluctuations, or scalar perturbations, is now well studied 
in the literature [l|, |2l , with a significant effort spent on 
improving the accuracy of analytical and numerical cal- 
culations (see, recent review in Q]). The non-Gaussian 
pattern of CMB anisotropics and polarization created by 
non-linear mapping associated with lensing angular de- 
flections aids the extraction of certain statistical proper- 
ties of the foreground mass distribution [^ . Weak lens- 
ing deflections by intervening mass also leak CMB polar- 
ization power in the E-mode to the B-mode |a]. This 
lensing B-mode signal presents a significant confusion 
when searching for primordial gravitational wave signa- 
tures in the CMB polarization [g. The lensing recon- 
struction techniques discussed in the literature, however, 
allow the possibility to "clean" CMB polarization maps 
and to search for a background of inflationary gravita- 
tional waves with an energy scale as low as 10^^ GeV 

0. 

Similar to gravitational lensing by density perturba- 
tions, if there is a background of gravitational waves in 
the foreground, then one would expect metric perturba- 
tions associated with these waves to distort and gravi- 
tationally lens background images [g. While the lens- 
ing deflections by the density fleld can be written as the 
gradient of the projected gravitational potential, lens- 
ing displacements due to gravitational waves can be de- 
cornposed to both a gradient and a curl-like component 
|3,II3jI11|- In these two components, gradient-type dis- 
placements are related to the lensing convergence, while 



curl-type displacements are related to the image rotation, 
though both types of displacements lead to image shear. 
While linear density perturbations do not produce rota- 
tions, second order corrections to weak lensing by scalar 
perturbations, such as due to the coupling of two lenses 
along the fine of sight, can produce rotational modes [ia |. 

While the study of CMB lensing by foreground density 
fluctuations is now well developed [3, the discussion of 
CMB lensing by foreground gravitational waves is lim- 
ited. In the context of large-scale structure weak lens- 
ing surveys with galaxy shapes [13 , the rotational power 
spectrum of background galaxy images when lensed by 
primordial gravitational waves in the foreground is dis- 
cussed in Ref. [lOj . In the context of lensing reconstruc- 
tion with CMB temperature and polarization maps, the 
curl component of the displacement fleld can be used to 
monitor systematics [ll| , though lensing by gravitational 
waves will leave a non-zero contribution to the curl com- 
ponent. 

Here, we extend the calculation in Ref. |lfl] and study 
both the curl- and the gradient-modes of the deflection 
field from primordial gravitational waves that intervene 
CMB photons propagating from the last scattering sur- 
face. Our calculations are both useful and important 
given the increasing interest on, and plans for, high sen- 
sitivity CMB anisotropy and polarization measurements, 
including a potential space-based mission after Planck, 
called CMBpol in the future. Such an experiment is ex- 
pected to study polarization B-modes in exquisite detail 
and it is important to understand potentially interest- 
ing secondary signals beyond those that are routinely 
mentioned in the literature. Based on the calculations 
presented here, unfortunately, we find that gravitational 
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FIG. 1: Lensing-deflection power spectra. Here, we show the 
gradient component from density perturbations (top curve), 
the curl (dot-dashed line) and gradient (dashed line) compo- 
nents from foreground inflationary gravitational waves. We 
have taken a background of gravitational waves with an am- 
plitude for the power spectrum corresponding to roughly a 
tensor-to-scalar ratio of 0.3 or a Hubble parameter during in- 
flation of 2 X 10" GeV. 



lensing of CMB by a background of primordial gravi- 
tational waves from inflation, with an amplitude below 
the current tensor-to-scalar ratio upper limit of 0.3, will 
produce an undetectable modification to anisotropy and 
polarization power spectra. Moreover, since the correc- 
tions are below the cosmic variance level, it is unlikely 
that one needs to account for these secondary corrections 
when making precise cosmological measurements. 

This paper is organized as follows. In Section m we 
discuss lensing by foreground gravitational waves by dis- 
cussing both the gradient and curl components of the 
displacement field. Section IIIII presents expressions for 
the weak lensing correction to the CMB anisotropy and 
polarization power spectra. We conclude with a discus- 
sion of our results in Section Hvl 



II. 



THE SPECTRUM OF EXPANSION AND 
ROTATION 



To establish the lensing correction to CMB anisotropy 
and polarization maps by foreground gravitational waves, 
wc first need to calculate the photon displacement on 
the spherical sky by gravitational waves in the fore- 
ground. We make use of synchronous coordinates and 
take the metric of a Friedman-Robertson- Walker cosmo- 
logical model as 



gp.i' 



-1 
I + H 



(1) 



where the scale factor is 0(77), and the conformal time 
is denoted by 77 with the value today of 770. Here, H is 



the transverse (V • H = 0), symmetric (H = H"""), and 
traceless (Tr H = 0) tensor metric perturbation associ- 
ated with gravitational waves, while I is the identity ma- 
trix. The photon propagation is governed by the geodesic 
equation 






(2) 



By changing variables to the conformal time 77 from the 
affine parameter A, the geodesic equation in the presence 
of foreground gravitational waves is [iU^ 



r-H-r r-(I + H)" 



.AH-iv.(r.H.r) 

. (3) 

where, to simplify the notation, we have not written out 
the explicit dependence of 77 in each of these terms. In 
here, the over-dot represents the derivative with respect 
to the conformal time. The full derivative d/drj can be 
separated to d/drj -I- r • V. Here, and throughout, S/h 
denotes the gradient that applies only to the metric per- 
turbation H; when not subscripted with 7f , the gradient 
should be interpreted as the one that applies to all terms, 
including the line of sight directional vector n. 

As gravitational fluctuation is very weak today, here- 
after we neglect H(?7o), and choose the initial conditions 
of the trajectory to be 



r(77o) = , r(77o) 



(4) 



we find that the general displacement on the celestial 
sphere induced by primordial gravitational waves is 



rvo 
r{ri) = n(77o - v) - H x 



H • n + 1(77' - 7?)[(n • H • fi)A - V H{n • H • n)] 



(5) 



(r;',x') 



where x' = (r/o — v')''^- Similar to Ref. [iQ, we have eval- 
uated the trajectory on the unperturbed path following 
the so-called Born approximation. One could potentially 
evaluate corrections to this approximation in terms of 
a perturbative correction to the path length, but these 
would be at the second order in metric perturbations and 
will be ignored. 

The total displacement can be separated to a part 
along the line of sight and a part perpendicular to it. 
The radial displacement leads to a time-delay effect, sim- 
ilar to the lensing time-delay associated with foreground 
potentials [iJI ■ This time-delay couples to the radial gra- 
dient of the CMB. due to finite extent of the last scat- 
tering surface, while the angular displacement couples 
to the angular gradient. As discussed in Ref. [I3, the 
overall correction to CMB anisotropy spectra from the 
time-delay effect is subdominant since the spatial gra- 
dient of the CMB is, at least, two orders of magnitude 
smaller compared to the angular gradient. There are also 
geometrical cancellations that make the time-delay effect 



smaller relative to angular deflections. Thus, we ignore 
the radial displacement and only consider the transverse 
component related to gravitational lensing angular de- 
flections. 

Using the transverse displacement, the angular deflec- 
tion projected on the spherical sky is A = [i"(??) — (n • 
r)n]/(77o — 77). These two-dimensional displacements can 
be related to usual quantities in gravitational lensing 
with the convergence and the rotation defined as Q 



Aa 



and 



_ {Aae''''):b 



(6) 



respectively. 

Similarly, we note that the general displacement on the 
celestial sphere can be decomposed to two components, 



^a = - T.C'fmyim-.a + hf^Y,^..,e\) , 



(7) 



hn 



where 



^'^ " l{l + V I ""'' 



^ ^dn>l™A„- = -^(7^ / ^ft^-'^ 






dhY,:Aa''e\ = T77^^ / dhY^l.LO, 



1(1 + 1) 



(8) 



where and (X) denote the gradient- and curl-type de- 
flections, respectively. For simplicity, we have dropped 
the dependence on the directional vector parameterized 
by n. The lensing by foreground density perturbations 
to the first-order only leads to a gradient-like displace- 
ment, while both components are generated when lensed 
by gravitational waves. We now calculate both the con- 
vergence and the rotational spectrum of the displacement 
field due to foreground stochastic gravitational waves. 



to two components and write K(n) = ki + K2 as 



«! = -;./ d,/p^^)(n.Vff)(n.H.n) 



2(^0 - Vs) 



/ dr]'{h • H • n). 



K2 



/,,, \Vo-VsJ 

1 rd^'[!?i^V^'_^^)(n.V«)2(n.H.n) 
4^^ \Vq-VsJ 



Va 



1 H" 

/ drj'i?/ - 77,)(n • VH)(n • H • n) , (10) 



where we have explicitly simplified the calculation by in- 
cluding terms associated with V • n. 

Note that we have also replaced rj ^ ris corresponding 
to the conformal time at the last scattering surface of 
CMB. To simplify, we decompose the metric perturbation 
into the Fourier component, 

1 r ^ 

H(x,^) = —-^ I d'ke''^-T{k,r,)Y,HjO,)e^ik){U) 

where we have introduced the gravitational wave transfer 
function that describes the time evolution of the metric 
perturbation with Hj(k, 77) = T{k,r])Hj{'k)ej{'k). 
The terms in equation H10() can be simplified as 



Kl 



^i/,(k)[n.e,(k).n](7k.n), 
i 

, 1 



■no-v 



rio-ri, 
3 



d-q 



2(770 - Vs) J,u (2^)'/' 

cT(fc,77')5]H,(k)[n.e,(k).n] 



d^ke*''-^' 



A. Gradient Spectrum 

Gradient deflections are associated with the expansion 
and, as defined above, can be described in terms of the 
convergence: ^(n) = — A'':a/2. In general the transverse 
divergence of a vector A can be rewritten as 

1 5 1 dA^ 

-(smMe) + 



sin 61 961 



smf 



V-I-(n-V)(n-I)-^^^ ) , (9) 



where r = rjQ — rj'. When substituting the form of r(77) 
from equation (5), the gradient terms here lead to terms 
that are due to V • n and V • H. We first consider the 
former and making use of the fact that difij ~ [Sij — 
fiinj)/{rjo — 77'), we separate contributions to convergence 



'' = L '' 4 (2.) 



''?' ''" 1 /"j3i„ikx' 



l^Jd^ke^^-'T{k,^') 



Vo - Vs 



^H,(k)[k2-(n.k)2][n.e,(k).n] 



1 /■'70 ^ I , 

dv'iv' ~Vs) .^ ...^ I rf^ke'*^-^ 



xT(fc,77') J2 H^m^ ■ e,(k) • n](zk • n). (12) 

j 

Here, x' = (770 — 77')n and e^ represents the symmet- 
ric, traceless polarization tensor that obeys Tr [ej(k) • 
efc(k)] = 2Sjk and k • ej(k) = with j in equation H12|l 
summing over the two linear polarization states. 

Since gravitational waves trace the wave equation with 



H - V^H + 2-H = mnGa^P 
a 



(13) 



where P is the tensor part of the anisotropic stress, say 
from neutrinos (see, Ref. [l^ for details); The term in 
the right hand side acts as a damping term for the evo- 
lution of gravitational waves and is important for modes 
that enter the horizon before matter-radiation equality, 
with a smaller correction for modes enter horizon after 
matter-radiation equality. Since these corrections arc 
not more than 30%, while the amplitude of the gravi- 
tational wave background is uncertain to more than or- 
ders of magnitude, we ignore such subtleties here as- 
suming no anisotropic stress; for a cosmological model 
dominated by matter, in Fourier space, one can write 
the evolution of H in terms of the transfer function as 
Tik,7j) = 3Mkrj)/{krj). 

We define the power spectrum of metric perturbations 
as 



(H.,ik)H*{k)) = (27r)3pr(fc)<5„<5(3)(k„k') 



(14) 



where we assume isotropy and equal density of gravita- 
tional waves in the two polarization states i and j. Fol- 
lowing Ref. |lCll | , we normalize the power spectrum to the 
Hubble parameter during inflation and take 






(15) 



Using this three-dimensional power spectrum for 
metric perturbations, the angular power spectrum of 
gradient-type deflections is 



C 



21 



m—l 

—1 E <l"-l 



(16) 



In equation. (16), we have also introduced an addi- 
tional correction to the gradient-type deflection spectrum 
with the term T®. As discussed in [13, this is associated 
with the gradient-type deflection pattern or convergence 
related to shearing of the last scattering surface by grav- 
itational waves present at that surface. To calculate this 
correction, we first evaluate the displacement vector at 
the last scattering surface 



H 



Hm - Vs) ■ 



(18) 



and then project this displacement vector to obtain the 
transverse displacement vector of 



(r„ • n)n 



Vo -Vs 



(19) 



This vector can be decomposed to the gradient- and curl- 
types deflections. Making use of the fact that the conver- 
gence is K{n) = —A'^.a/2 and taking the Fourier trans- 
forms, we write the required term in the equation. (16) 
as 

k 

T® = -i^{r]o-Tls)T{k,ri,)[d^{x-'^ji{x)%^k{ri^_^,)] 

-T{k,ris){x~'^ ji{x)%^k(na-v.) ■ (20) 

This completes the calculation of gradient-type deflection 
power spectrum by taking into account the lensing by in- 
tervening gravitational waves between the last scattering 
surface and the observer and the metric-shear correction 
related to gravitational waves present at the last scatter- 
ing surface from which photons propagate. 



7n—l 



7n—~l 



dhYC^{nMh)\') 



(1 + 2)1 



PH + l)"^ {1-2)1 
where the terms are again 



f d^kPTik)\Ti+T2+T; 



2 
m I ' 



Ti 






T(fc,7?') 



X [d^{x '^jl{x))\^=kir,o-ri') 

3 



To 



VO - Vs Jr,, 

'2 



d7fT{k,jf){x ^ji(a;))U^fc(,)o -»)')' 



dv'(^^^^^]{v'-Vs)T{k,r,') 



'no -Vs. 

X [{^+dl){x-'^]l{x))\,=k{r„^r,')\ 

2k 



drj'{rj'-ijs)T{k,if) 
Vo - Vs J,,, 

X [dxix~^jlix))U=k{r,o-v')] ■ 



(17) 



Note that the convergence power spectrum due to gravi- 
tational wave deflections is Cf' = P{1 + 1)^(7/'® /4. 



B. Curl spectrum 



Rotation is defined to be w(n) = ^{Aae'^^);b, which 



leads to 



w(n) = --n- (V X r(n,77s)). 



(21) 



Equation © gives 



U! 



1 r° 

- j dr^' [n • (V X H) • n] , (22) 



since V x n =0. Here we define the curl of the second 
rank tensor H by (V x H)^; = e^jkdjHki- 

The curl-type deflection spectrum from this term is 



cf 



m— / 

m^ — l 

m—l „ 

1)12(1 + 1)2 E i\ J dnYi:nif^)uj{n) 

7n— — l 



{2l + l)P{l + l) 

TT {l + 2)\ 

P{l + l)^{l-2)\ 



ld^kPTik)\Ts + T^\',(23) 



where 



T. 



2k 



diT{k,7^'){x-^ji{x))U^k^^^_^') , (24) 



Again, one can write the power spectrum of rotation as 
C^^ = P{l + lYcf /A. 

Here also we include the correction to the rotational 
spectrum due to metric perturbations at the last scat- 
tering surface. Similar to convergence, by following the 
same procedure as before but taking l/2([A„i]ae°''):(, of 
Equation H19|) . we get 



rp^ 



HVO~Vs)T{k,r]s){x '^jlix))\a:=k(rjo-Vs)- (25) 



For comparison, we note that density fluctuations 
along the line of sight lead to gradient-type deflections 
only. The resulting contributions are described in terms 
of the angular power spectrum of projected potential Cf 
that is well studied in the literature Q . We do not repeat 
those derivations here, but will provide a comparison of 
lensing under gravitational waves and lensing by mass in 
the discussion later. 



III. LENSING OF CMB BY GRAVITATIONAL 
WAVES 

In this Section, we will discuss the analytical calcula- 
tion related to how foreground gravitational waves mod- 
ify the background CMB temperature anisotropy and po- 
larization patterns. Here, we will concentrate on the an- 
gular power spectra of CMB observables. The lensing of 
CMB by foreground density fluctuations is formulated in 
Ref. I2I and we follow the same procedure here. Since de- 
flection field power spectra peak at large angular scales, 
we present an analytical formulation as appropriate for 



the spherical sky. The calculation we consider here, how- 
ever, is perturbative and one expects important correc- 
tions beyond the first order in the deflection angle. For 
foreground gravitational waves, such issues can be ig- 
nored as the overall modification to the anisotropy and 
polarization spectra is small. 



A. Temperature Anisotropies 

Following Ref. (2|, the lensed temperature field 9 can 
be expressed as 



(n) = e{h + A) = e{n) + V^e* • A„ + -V''V"6i • AaAb 



(26) 



where 0{n) is the unlensed temperature fluctuation in 
the direction n. The temperature field can be expanded 
to multipole moments such that ^(n) = J2im^hnYim{i>)- 
Taking the spherical harmonic moment of equation H2t)|) 
and using equation ((7J, we find 



9i^ + / dhY, 



V"0 • A, + -V''V'^0 • A„Ab 



q _ Y^ / 7-© /) I,© 

Iimil2m2 



+1, 



Imlimil2m2 '■i"^! l2m.2 



- V ( I® e, h® h- 

2 Z^ y'lmlimil2m2hm3^limill'i^„,^ll'i 



2m2 lams 



Iiniil2m2hm3 



^■^ Imlimil2m2l3m3^ hmii '■l2m2 H^mi j i 

where the integrals are 



(27) 



Irnlimil2rrt2 






f rinY* y-" Yi T^ 

•^Imhmihmzhma ~ J "'^Mm^iiTni^'2m2:aM3m3:& ' Jlmhmil2m2l3ma ~ / "■^■Mm^iimi^2m2:cJj3m3:d^ a^ b- (28) 



The lensed temperature anisotropy power spectrum spectrum is 

Ct ^CUY: Cl (Cf S® + Qf Sf) + Cf Y: « S® + Qf Sf) , 



ii;2 



(29) 



where 



s® = 


- V 1/® P 




mim2 


sf = 


/ ^ \ Inilimil2m2^ 




mim2 


^2® = 


~ 2 / .ImlmUrmUrm + C-^- 




mi 


sf = 





and c.c. is the complex conjugate. The terms 5® and 
S® are similar to those involving lensing by foreground 
(30) density perturbations Q- First, the integral I^i ^ ; ^ 
can be simplified through integration by parts and noting 
V^yjm — —i{l + i)Yi„i and the general integral of three 



10 





FIG. 2: The lensing modification to CMB power spectra for density perturbations and for gravitational waves. Left: Tempera- 
ture fluctuations. The top curve is the primordial power spectrum. The middle curve is the secondary anisotropy contribution, 
\Ct — C;| , to the temperature power spectrum from lensing by density perturbations, and the lower curve is for lensing by 
gravitational waves, assuming the maximum IGW background consistent with current data with the deflection power spectra 
shown in Figure 1. Right: Temperature-E polarization cross correlation, with curves following the left panel. In both panels, 
thin long-dashed line is Ci/l at each multipole with Ci related to the intrinsic anisotropy spectrum; while the cosmic variance 
is Ci/Vl, Ci/l denotes the level at which one must control the systematics, if the effects resulting systematics apply to a wide 
range of multipoles. The lensing by density fluctuations cannot be ignored as the corrections are well above the cosmic variance 
limit and will be detectable in upcoming anisotropy data. The lensing by foreground gravitational waves, however, are below 
the cosmic variance limit, suggesting that they will remain undetectable, but above the systematic level when I > 10'^. 



spin-spherical harmonics over the sky: 






{2h + l){2l2 + l){2h + l) 



h 

Sl -S2 -S3 



47r 
W h h 

-7711 "72 TO3 



(31) 



when Sl = S2 + S3. We note that under parity inversion, 
sYim -^ {—^Y-sYim,, which is a useful property when we 
discuss lensing modifications to the CMB polarization 
field. 

With Si = and noting that o^Zm = ^Zmi 



/; 



iml.m.l.m. = ^^Cl + 1) + Uh + 1) - 1(1 + 1)](-1)^ 



(2; + l)(2;i + l)(2Z2 + l) 



/ h h 



47r 

I h h 
—m mi TO2 



0'* — "^ "^- ■^" ' ' ^ ' 

Using the orthonormality relation of Wigner-3j symbols 



E 



mini2 

we can write 



h h h 
TTii 7712 m^ 



h h h 

Till 7712 1^3 



1 



2^3 + 1 



5® = 



21 



— fF® V 



(33) 



(34) 



with 



pa 



^[h{h + l)+l2{l2 + l)-l{l + l)] 



^ J {2l + l)(2h + l){2l2 + l) fl^ h h , ^ (35) 



All 







Though tedious, the calculation related to Sf can be 
simplified using the gradient relation for s phe rical har- 
monics by raising and lowering of the spin |l7l | : 



^Yim = \r-^^-^[iYiram+ - _ii^™m_], (36) 



where 



m± = ^(eeT«e^)- 



(37) 



Combining these derivatives with the general integral in 
equation 131|) leads to 



Iml 



,™,z.™. = -T^^h{h + i)Hk + i){-iy 



{2l + l){2li + l){2l2 + l) 



An 
i ti '2 \ / ^ ^1 ^ 



-771 Till 7712/ \0 



\ ^A [1 - (-l)'+'^+'^](38) 



Again using the orthonormality relation in equa- 
tion (123), 



^ ^21 + 1 ^ "1'^^ ' 



(39) 
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FIG. 3: The lensing modification to CMB polarization spectra 
for density perturbations and for gravitational waves. From 
top to bottom, these curves are for primordial power spec- 
trum in the E-mode (solid line), lensing correction to E-mode 
from density fluctuations (long -dashed line), lensing correc- 
tion to B-mode from density fluctuations (short-dashed line 
labeled "Scalars"), lensing correction to B-mode from gravi- 
tational waves (short-dashed line labeled "Tensors"), and the 
lensing correction to E-mode from gravitational waves (long- 
dashed line). For reference, the thin dot-dashed line is the 
systematic level of the B-mode lensing power spectrum, Cijl, 
from the density fleld. While the corrections from lensing by 
foreground gravitational waves is below the cosmic variance 
limit of the lensing B-mode power spectrum and, again, un- 
detectable in anisotropy maps, they may become a source of 
systematic in all-sky maps with no instrumental noise and 
other secondary signals and foregrounds. 



with 



\ = lVhih + iMi2 + i)(l^ !\ \') (40) 

[!-(-!)' 



{2l + l)i2h + l){2h + l) r, , ,,,+,,+,, 



Att 



To calculate Sf and Sf, we first note that 



/ , ^'imi:a-'/imi:h 



^i(/i -H)(2^i + 1) 



where gab is the usual metric of unit sphere, 



_ 1 
^"^^ ~ ' sin^ 



9ab (41) 



(42) 



These allow us to show that 

2^1 + 1 



mi 

4{l + l)h{h + l)- 



Stt 



(43) 



Also 5*2 = S®- Finally, combining all expressions, we 
can write 

Cf ^ Cf-lil + l)RCf 



E 



'1 

21 + 1 



C«J^ + C«J^ .(44) 



where F®j and Fjf ^ are given in equations (|35|1 and 
(40), respectively, and 



^ = E^i(^i + l 



2?i + 1 r 






(45) 



This expression is similar to that of equation (62) of 
Ref. [31 when C/* = and C/* is identified as the power 
spectrum of projected lensing potentials due to interven- 
ing density perturbations between us and the CMB. 



B. Polarization 

The lensing effect on CMB polarization can be de- 
scribed similar to temperature anisotropics by mak- 
ing use of the remapping ±X{h) ~ ±-^(n) + 
Vi(/)(n)V*±X(n) + ...+, where ±X = Q ± iU, where 
we have simplified the notation by replacing the spin- 
dependent gradients with a covariant derivative that acts 
on the spin-components of the symmetric tensors that 
are traceless. Here, we have also ignored the rotation 
needed to align the polarization basis vectors between 
the lensed and unlensed fields, but this rotation is unim- 
portant when considering displacements along the lines 
of constant azimuthal angles jl8j . We refer the reader to 
Rcfs. [isl [13 for details of our shorthand notation and 
why it can be used for the lensing of the polarization 
pattern on the spherical sky. While our notation here is 
simple and follows that of Ref. [J , the final result is the 
same from what one gets by using the standard notation 
of differential geometry. This is due to the fact that the 
overlapping integrals involving spin harmonics that we 
will perform remain consistent with our simplified nota- 
tion. 

As is well known, the CMB polarization components 
form a spin-2 field and are expanded in terms of the 
spin- weighted spherical harmonics such that ±X{n) = 
X)zm ±^i'n±2^m(A). Instead of Stokes parameters, the 
more popular, E and B modes are are given by ±Xim = 
Ei„i ± iBim.. We will discuss lensing modifications to 
angular power spectra of E- and B-modes as well as the 
cross-correlation between E and 6, while cross-correlation 
between E and B, and between B and 9 are ignored as 
these are zero through parity arguments. Furthermore, 
the modifications to polarization by lensing do not vio- 
late parity conservation. 

Taking the spherical harmonic moment of the polar- 



ization field under lensing, we write 



±Xlm — ±Xlm — 

E 



Iim\l2m2 






After straightforward but tedious algebra, the lensed 
power spectra of iJ-modes, i3-modes, and the cross- 
(46) correlation between iS-modes and 9 are 



+±2/; 



hf. 



E 



Iral\7n\l2ra2 l2'n^2 



±^lirni \±2'J Imlimil2m2l3m3"'l2m2"'hmri 



1 70 iM l8* 

=t2 I'mlimil2m2l3ni3 l2m2 l^ms 



where 



rffi 
^2 hnlim,il2 



= j dh{±2YC„,){±2Y(^^jYi,m2:a 
±2Ji,nlimil2Vi2l3m3 ~ / ^"^(±2^/™) (±2^iimi )^'2m2:a5^i3,„3:6 

±2lZihrail2rn2 ^ / ^"(±25^^) (±2>"iimi )^i2m2:&e''a 






(47) 



/ dn(±2y4) (±2y;i™jYi2m2:cy;*m3:de''ae''b. 



c,^ 



cf 



^SB 



Q'^ + ^ E [<^Sf + cf.Sf] [{Ct + CE,) + (-l)^(Qf - Qf )] + cPY.(<^Sf + cf.Sf), 
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-\Y.(^ + (-l)^)(cf 02^? + <02^f )Qf + \cr E (< (2^2® + Sf) + Cj^: i^Sf + ^f ))(48) 



cr = c'l 



where L = I + li + I2 and 

2'S'f = 2^ l±2-^®i;imi 



i2"l2l ' 



(49) 



c® _ Y^ I , /""^ |2 

2-^1 — 2^ l±2-'imiimii2m2l ' 
rnirn2 

C© _ 1 V^ , T-e I „ „ 

mi 

■^2 ~ o / ,±2'JimlmUm,hrm "t" CX., 

T7ll 

c© _ y^ (T® J® \ 

02^1 / J \ Iml-ira-il2inri2 +2 Iml-im-il2m2^ ^ 

77117712 

o® _ \^ (1-0$ r® ^ 

02'-'l / y \ Im,limil2m2 +2 /m/i7nii27n2'' 

7^17712 

To simplify terms in equation (49), we again make 
use of the integral relations outlined earlier when de- 
scribing lensing of temperature anisotropics. In the 



r 



case of polarization, these relations need to be general- 
ized for integrals over spin-weighted spherical harmon- 
ics. First, the integral related to the gradient spec- 
tra is straightforward. Making use of the fact that 
V^±2Y/7n = [-^(^ + l) + 4]±2^m and usiug cquatiou (|3lJ, 
we find 

C© _ ^ I p© |2 



^^^^I i2l + l)i2h + l){2h + l) fl^ I, l,^ I ^^^^ 



Att 



2 0-2 



This is exactly the relation that one encounters when 
lensing the polarization field by foreground density fluc- 
tuations |3|- 

The integral related to the curl-type displacement is 
tedious, but can be simplified using relations involving 
raising and lowering of the spin and gradient of the spin- 



weighted spherical harmonic. To calculate 



2 Iral\7n\l2rri2 

we note 

m_ • Vsllm 

and the relation 



fdni±2Y;j±2Yl,^,..aYi,rn._:be''^51) 



{l-s)il + s + l) 



-s+i-i;™ 



i1. 



{l + s){l-s + l) 



(m+),(m_)j + (m_),(m+)j == gy, (53) 

to write V sYim = (m+m_ + m_m+) • WgYim as 



.Yr 



il^s){l + s + l) 



s+1 



Yhnin^ 



{l + s){l-s + l) 



-lYhnm^. 



(54) 



This leads to 



■^hmi 



irathm-i ~ / C?n(+2i^m)+2^imi :a^2m2 ib*^ " 



-iFz™ , (52) = z(-l)' 



^2(^2 + 1) /(2;+i)(2Zi + i)(2;2 + i) 



47r 



-m TTil 7Ti2 



(/ 



(/i + 2)(;i-l) (I I, l2 
2 -1 -1 



{h-2){h+3) (I W h 
2 -3 1 



such that 



(55) 



c«> _ Y^ I r® 

2"^! — /_^ I'^^lmhmil 



1 



2/ + l" "i'= 



2-^!ii/2 



2 -1 -1 



r 



Furthermore, with s = ±2, equation H43() can be gener- perature, we find 
alized to 



(56) 



Z2(^2 + l)(2/ + l)(2;i + l)(2?2 + l) / /(?i+2)(/i-l) /; li ?2\ [h - 2){li + 'i) I li I2 
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02 



Qtg> 
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21 



-jiPihi2)i+'^^ihhy 



(59) 



and as in the case of lensed temperature anisotropics, 

+2S2 = +2S2 ■ (58) 

Putting the terms together, we can write the lensed 
For the cross-correlation between E-modes and the tem- power spectra in polarization as 

I 



Q" 



c^- 



f -if + l~ mcf + ^^^ E [ct'i^Fl.f + Cf (2< J^] [(Qf + Qf ) + (-f)^(C,f - Qf )] 

'l'2 

Vtt E K(2< J' + Qf (2< j^] [(Q? + Qf ) - i-incF^ - CO] 



Cf = Cf - (/2 + ; - 4)i?Cf + 



2(2^ + 1) 



Q 



BE _ r<SE 



cr -{I' + l 






ci 



SB 

^2 ' 



(60) 



with R from equation (|45|l . 



The case of CMB Icnsing by foreground density fluctu- 
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ations is simply the replacement of C/* with the power 

spectrmn of projected potentials and with Cj^ — 0. The 
lensing contribution to the B-niode from the foreground 
density field [5J act as the main contaminant in detect- 
ing the primary gravitational wave signal in B-modes of 
polarization. To see the extent to which lensing by grav- 
itational waves themselves may become important in B- 
mode polarization studies, we will assume CP to be zero 
and present a comparison between CP from density fluc- 
tuations and CP from gravitational waves. As we find, 
the secondary lensing from gravitational waves is smaller 
than the cosmic variance level of lensing B-modes from 
density perturbations and will remain undetectable in 
anisotropy maps. While below the cosmic variance limit, 
Ci/Vl, for anisotropy measurements, systematics must 
generally be controlled to a level far below this; if sys- 
tematics apply to a wide range of multipoles, then one 
must control their effects to Ci/l. For lensing by fore- 
ground gravitational waves at the maximum amplitude, 
we find that the corrections are above this level when 
/ > 10"^ suggesting that one only needs to be concerned 
of these signals in all sky maps with no instrumental noise 
and other secondary signals. 



IV. RESULTS AND DISCUSSION 

In Figure 1, we show a comparison of C/* , Cj^ , and 
the angular power spectrum of deflection angle from pro- 
jected density perturbations along the line of sight to 
77s at a redshift of 1100 corresponding to the CMB last 
scattering surface. In calculating the power spectra of 
lensing from foreground gravitational waves, we have as- 
sumed an amplitude for the tensor modes with a value for 
Hj in equation (|15|l of 2 x 10^^ GcV. This corresponds to 
a tensor-to-scalar ratio of 0.3, which is roughly the upper 
limit allowed by current CMB and large-scale structure 
observations |20l| . 

The curl spectrum of deflections from gravitational 
waves has been previously discussed in the literature in 
the context of weak lensing surveys with galaxy shapes 
[lOj. The gradient- type displacement spectrum from 
gravitational waves discussed here is also important and 
cannot be ignored when calculating modifications to 
CMB temperature and polarization anisotropics. Note 
that Cj^ and Cj' peak at large angular scales corre- 
sponding to i! = 2 to ^ = 10. To compare lensing by 
gravitational waves and lensing by mass, we calculate 
the rms deflection angle through 0^,jjg — [J2 l{^ + 1)(2' + 
l)/47r(C/' +C'/' )]■ For gravitational waves related spec- 
tra shown in Figure 1, 0rms is 3.62 x 10~^ radians or 
roughly 7 arcsecs, which is a factor of 20 smaller than 
the rms deflection angle for CMB photons under density 
fluctuations, where ^rms ~ 7x lO"**. The coherence scale, 
where the rms drops to half of its peak value, is about 
'^ 60 degrees for lensing by gravitational waves, while 
for density perturbations the coherence scale is about a 



degree. With such a large coherence scale and a small 
rms deflection angle, foreground gravitational waves de- 
flect large patches of the CMB sky by the same angle of 
about 6 arcsecs, resulting in an overall small modification 
to the CMB anisotropy and polarization spectra, when 
compared to the case with density fluctuations alone. 

The differences between lensing by density perturba- 
tions and lensing by gravitational waves is clear in Fig- 
ure 2, where we show modifications to the tempera- 
ture power spectrum and the cross power spectrum be- 
tween temperature and E-modcs of polarization. The 
secondary lensing correction from the foreground gravi- 
tational waves is smaller than the cosmic variance level 
of intrinsic CMB anisotropics. Even with perfect CMB 
observations devoid of instrumental noise, it is unlikely 
that the lensing modification by primordial gravitational 
waves in the foreground of CMB will be detectable. As 
shown in Figure 2, however, the corrections are above the 
systematic level of Ci/l for primordial anisotropy mea- 
surements when / > 10"^. Thus, in the extreme case 
where one is dealing with perfect all-sky maps cleaned 
of foregrounds and other secondary signals, one could be 
concerned that these effects act as a source of system- 
atic error for primordial anisotropy measurements. As is 
clear from Figure 2, while lensing by gravitational waves 
is not significant, lensing of the CMB by foreground den- 
sity perturbations will be detectable as the modifications 
are well above the cosmic variance limit. 

In Figure 3, we summarize our results related to lensing 
of polarization anisotropics in terms of the power spectra 
or E- and B-modes. Again, the lensing effect by fore- 
ground gravitational waves is below the cosmic variance 
level of the dominant signal in the E- and B-mode maps, 
but above the systematic level. In the case of B-modes, 
lensing by density perturbations will remain the main 
contaminant in searching for gravitational wave signa- 
tures from the primary B-mode power spectrum. Lensing 
signals by foreground gravitational waves are unlikely to 
affect reconstruction techniques that attempt to remove 
the lensed B-modes when searching for a low amplitude 
gravitational wave background Q . 

While we have only considered angular displacements 
on the sky, gravitational waves also contribute to a vari- 
ation along the line of sight that can be described as 
a time-delay effect. Just as angular displacement cou- 
ples to the angular gradient of the CMB, the radial dis- 
placement couples to the radial gradient of the CMB. 
These effects, however, are smaller due to lack of radial 
structure in the perturbations that form the primordial 
anisotropy spectrum in the CMB, such as the acoustic 
peaks. There are also geometric cancellations associated 
with the projection of line-of-sight time-delay modula- 
tions to an anisotropy pattern on the CMB sky |15l| . 
Thus, it is unlikely that our conclusions related to lensing 
by foreground gravitational waves are affected by includ- 
ing the time-delay effect. 

In general, our results are consistent with those of 
Ref. [13 who studied the possibility of measuring the 
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gravitational wave background amplitude using weak 
lensing surveys of galaxy shapes and using the curl-mode 
of the shear. Even with an optimistic survey with a large 
surface density of galaxies to measure shapes, the grav- 
itational wave signal in the shear remains undetectable 
below the noise. 

To summarize our calculation, while weak lensing dis- 
tortion of the cosmic microwave background (CMB) tem- 
perature and polarization patterns by foreground density 
fluctuations is well studied in the literature, we noted the 
lack of a detailed description related to lensing modifica- 
tions by foreground gravitational waves or tensor pertur- 
bations. Here, we have presented an analytical formula- 
tion on how CMB anisotropics and polarization patterns 
are distorted by a stochastic background of primordial 
gravitational waves between us and the last scattering 
surface. Our analytical formulation is useful when study- 
ing general lensing of any background source by fore- 
ground gravitational waves. 

While density fluctuations perturb CMB photons via 
gradient-type displacements only, gravitational waves 
distort CMB anisotropics via both the gradient- and the 
curl-type displacements. The latter can be described 
as a rotation of background images in the presence of 



foreground gravitational waves while the former is re- 
lated to the lensing convergence. For a primordial back- 
ground of gravitational waves from inflation with an am- 
plitude corresponding to a tensor-to-scalar ratio below 
the current upper limit of ~ 0.3, the resulting modifica- 
tions to the angular power spectra of CMB temperature 
anisotropy and polarization are below the cosmic vari- 
ance limit, but above the systematic level. Thus, it is 
unlikely that planned high sensitivity CMB observations 
warrant an accounting of the secondary contributions dis- 
cussed here as they are not expected to affect precise pa- 
rameter measurements; if observations are all-sky mea- 
surements with no instrumental noise, then these effects 
may be present in the form of systematic corrections to 
the primary anisotropy and polarization measurements. 
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